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Coherent state functional integral for the minisuperspace model of loop quantum cosmology is 
studied. By the well-established canonical theory, the transition amplitude in the path integral 
representation of loop quantum cosmology with alternative dynamics can be formulated through 
group averaging. The effective action and Hamiltonian with higher-order quantum corrections are 
thus obtained. It turns out that for a non-symmetric Hamiltonian constraint operator, the Moyal 
(star) -product emerges naturally in the effective Hamiltonian. For the corresponding symmetric 
Hamiltonian operator, the resulted effective theory implies a possible quantum cosmological 
effect in large scale limit in the alternative dynamical scenario, which coincides with the result 
in canonical approach. Moreover, the first-order modified Friedmann equation still contains the 
^ ' particular information of alternative dynamics and hence admits the possible phenomenological 

i distinction between the different proposals of quantum dynamics. 
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I. INTRODUCTION 

In the last two decades, considerable progress has been made in loop quantum gravity (LQG), which is a background 
independent approach to quantum gravity [1-4] . The starting point of LQG is the Hamiltonian connection dynamics 
of GR rather than the ADM formalism. By taking the holonomy of sw(2)-connection A l a and flux of densitized 



triad E b , as basic variables, the quantum kinematical framework of LQG has been rigorously constructed, and the 



• Hamiltonian constraint operator can also be well defined to represent quantum dynamics. Moreover, a few physically 
i-H \ significant results have been obtained in the minisuperspace models of loop quantum cosmology (LQC) [5, 6]. The 
i-H . most interesting one is the resolution of big bang singularity in LQC [7-10]. Besides the canonical formalism, the 
' so-called spin foam models were proposed as the path integral formalism of LQG [1]. However, whether the two 
approaches are equivalent to each other is a longstanding open question. Thanks to the development of LQC, we 
have a much simpler theory to address this question. As symmetry-reduced models, there are only finite number of 
degrees of freedom in LQC. Following the conventional method in quantum mechanics, one can find the path integral 
formalism of LQC starting with the canonical formalism. This approach has been implemented by a series of works 
[11-14] with the scheme of simplified LQC [15]. Here one employed the complete basis of eigen-states of the volume 
operator to formulate a path integral with somehow descrete-steps, which inherited certain properties of spin foams 
[16]. Moreover, the first-order effective action for the path integral was derived by this approach [13, 14], which 
implied the origin of singularity resolution of LQC in the path integral representation. In canonical LQC, the effective 
Hamiltonian constraint with higher-order quantum corrections could even be obtained by the semiclassical analysis 
using coherent states, which implied a possible effect of quantum gravity on large scale cosmology [17-19]. It is thus 
interesting to see whether the effective Hamiltonian can be confirmed by some path integral representation. Since the 
higher-order corrections of the Hamiltonian come from the quantum fluctuations, a natural attempt to achieve them 
is to employ coherent state path integral [20]. 

In LQC, the Hamiltonian constraint equation is usually presented to Klein-Gordon like equation by coupling with 
a massless scalar field, where the corresponding gravitational Hamiltonian operator, as some multiplication of several 
self-adjoint operators, is non-symmetric in the kinematical Hilbert space [9, 10, 15]. While this treatment is essential 
in order to obtain the physical states satisfying the constraint equation, it also provides elegant physcial models to 
examine the so-called Moyal ^-product in quantum mechanics. At the very beginning, Moyal proposed the ^-product 
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in order to clarify the role of statistical concepts in quantum mechanics system [21]. Then this idea were generalized 
to many situations including quantum spacctime itself. In canonical quantum theories, the * -product can also be 
understood by coherent state approach [22]. Thus it is also possible and desirable to derive the ^-product in coherent 
state functional integral approach. This idea has been accomplished in the WDW quantum cosmology and LQC [23]. 
In these models, the Euclidean and Lorentz terms in the gravitational Hamiltonian constraint are combined together 
since they are proportional to each other in spatially flat and homogeneous cases. However, this is impossible in the 
full theory, where the Lorentz term has to be quantized in a form quite different from the Euclidean one [3, 4]. This 
kind of quantization procedure which kept the distinction of the two terms was proposed as alternative dynamics 
for LQC [18]. Hence we will study the coherent state functional integral in spatially fiat isotropic FRW cosmology 
coupled with a massless scalar held (j> in the alternative quantization framework. 



II. COHERENT FUNCTIONAL INTEGRALS 



We consider the following Hilbert-Einstein action of gravity coupled with a massless scalar field: 

S = J d i x^R+ l - J #xy/=g<r<t>, v <t>,v (2.1) 

In the spatially flat model of FRW cosmology, we fix a space-like sub-manifold S, which is topologically R 3 and 
equipped with Cartesian coordinates x l {i — 1,2,3) and a fiducial fiat metric °q a b- The physical 3-metric q a b is then 
determined by a scale factor a satisfying q a b = a 2o q a b- It is convenient to introduce an elementary cell V and restrict 
all integrations to this cell in Hamiltonian analysis. The volume of V with respect to °q a b is denoted as V a . As in 
the full loop quantum gravity, we employ the new canonical variables (A l a ,Ef) [24]. Due to the homogeneity and 
isotropy, we can fix a set of orthonormal cotriad and triad , °e") compatible with °q a i, and adapted to V. Then 

— (1/3) — (2/3) 

the basic canonical variables take the simple form A l a = cV uj l a , Ef = p^/°qV ef and thus are reduced to 
(c,p) with the Poisson bracket: {c,p} = 8irG-f/3, where 7 is the Barbero-Immirzi parameter. Following the fl-schem 
of "improved dynamics" [10], the regulator £l used in holonomies is given by fx = \J ~K]\p\, where A = 4\/3Tr"/£ 2 is 
a minimum nonzero eigenvalue of the area operator [6]. In order to do the semiclassical analysis, it is convenient to 
introduce new dimensionless conjugate variables [15, 18]: 

flc sgn(p)|p|§ 

V ' = ^7Z> (2 ' 2) 

with the Poisson bracket {v,b} — — j- , where the Planck length £ p is given by £ 2 = Gh. From the matter part of 

action (2.1), we can get p<f, — a v ^ and the poisson bracket: {(j),p^\ = 1. The kinematical Hilbert space of the 
quantum theory is supposed to be a tensor product of the gravitational and matter parts. In LQC, one employed the 
standard Schrodinger representation for matter to construct Hilbert space WJ-jf * , while gravity was quantized by the 
polymer-like representation [24]. Thus quantum states in the gravitational Hilbert space of LQC are functions ^(v) 
of v with support on a countable number of points and with finite norm \\^\\ 2 := ^2 V \^(v)\ 2 [25]. Hence the inner 
product is defined by a Kronecker delta (v'\v) = 5 v > tV . The basic operators act on a quantum state ty(v,<f>) in the 
kinematical Hilbert space "Hf[„ v as: 

vty(v,<f>) = v^{v,(f>), e lb ^{v,4>) = ^(v + 1,0). (2.3) 
To obtain the physical states, one has to solve the quantum Hamiltonian constraint equation: 

-d-*(v,(t>)=^-j± + ej9(v,<l>) = 0, (2.4) 

where & = 8e + ©l is a positive second-order difference operator defined by: 

9 E = [v(v + 2)fr(u + 4, (/)) - 2v 2 ^{v, <P) + v(v - 2)V(v - 4, 0)] , (2.5a) 

§ L • *(w, 4>) = ~ 3?rG ^ 6 +7 - i v ( v + 4 )*( w + 8 ^) - 2 « 2 *( u , 4>) + v(v - 4)*(w - 8, 4>)} . (2.5b) 

Here we use the alternative quantization scheme proposed in [18] in which the Euclidean and Lorentz terms in the 
gravitational Hamiltonian constraint are treated separately. Together with the simplified treatment in [15], we can get 
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the operators in Eqs. (2.5a) (2.5b) corresponding to the Euclidean and Lorentz terms respectively. Solutions to the 
constraint equations and their physical inner products can be obtained through the group averaging procedure. Now 
we concern about coherent state functional integrals. The (generalized) coherent state of the matter part is labeled 
by a complex variable z a :— —^{fto + f " 2 ? 3 ^) an d defined by 



* Zo ) := / d<P e-^^'e^'*-*"^), (2.6) 



which is the eigenstate of the annihilation operator z — -ji^{4> + jiP<t> ~ 2 )i where a describes the width of the wave- 
packet or quantum fluctuation. It satisfies the key properties of a coherent state, namely, saturation of Heisenberg's 
uncertainty relation, resolution of identity and peakness property. On the other hand, due to the polymer-like 
structure, the coherent state of LQC is different from that of the matter part. Here one can define ( D = -^( v o + ib d 2 ) 
to label the generalized coherent state [17, 25]: 



(* Co \:=J2e~ i ^e- ib ° {v - V ° ) (vl (2.7) 



where d is the characteristic width of the wave packet and 1 <C d <C v Q because of the semiclassical feature. For 
practical use, one defines the projection of this state on some lattice of variable v, saying the shadow state [25]: 



^ Co ) shad : = e- {± ^e ib °( k ~ v ° ) \k), k e Z, (2.8) 

k— — oo 



where we chose the regular lattice {v — k, k <G Z}. This shadow state also has the analogous properties of a coherent 
state. The resolution of identity now reads 

where the identity I is in the subspace in which the states have support only on the regular lattice. The whole coherent 
state of LQC reads |* z J|* Co ) = |* Z J <g> |* Co ). 

In the path integral of the conventional non-relativistic quantum mechanics, one needs to compute the matrix 
element of the evolution operator e ~ lAtH within the time interval At. However, the situation of cosmology in GR 
is very different, since we are considering totally constrained systems and the operator C is not a true Hamiltonian. 
Instead, we start from the physical inner product, i.e., the transition amplitude, of coherent states with normalization: 

A([*f],[Vi]) = Vl ,,,,,,,, „' — . (2.10) 

U /M U II^IIII^IIII^JHI^II 1 ' 

To calculate the transition amplitude, we split a fictitious time interval At = 1 into N pieces e = . To deal with the 
parameter a in group averaging procedure, we employ the trick in [14] to generalize the one single group averaging 
to multiple ones, i.e., 



lim / 



a n 



dae laC \^ M n) (2.11) 

lim — - / da N --- ■— \ da 2 da ie ^ + - +a ^ c \^ kin ), V |* kin ) eMMn. (2.12) 

■ ,oci ^oo 2a No J _ & 2a 2o J^ &2n J- & . 



In order to trace the power for expansion, we re-scale the parameters by a n — ea n (n = 1, • • • ,N) and thus rewrite 
the exponential operator as: e J ^-™=i tanC = Y[ n =i e ltanC . Inserting N times of coherent states resolution of identity 
of |* Zo ) and Eq. (2.9), Eq. (2.10) can be castcd into 

A([9f], [**]) = lim -J- I*" da N ■ ■ ■ J- C° da 2 ■ e C° d ai A% a \ (2.13) 
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where 
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„ 27rft "' 2tt?i 

00 n=l 
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l* 2 



dvN-i ■ ■ ■ dv\ 



71 db N -! dh -A- (*„„|e 



2tt 



2tt 



n 

n=l 



l* 2 



Zn-1 
-iea n 



1*^1111*^-11 



(2.14a) 



(2.14b) 



with ztv = , zo = z i 1 W = Vf > an( i Vo = Vi- Notice that the characteristic widths cr and <i at different steps are not 
necessarily the same. So we have to denote a n and d„ in the semiclassical states |* Zn ) and \^Q n ) respectively at the 
"n-step" . Now the main task is to calculate the matrix elements of the exponential operators on coherent states. The 
exponential operator e l£Q " c can be expanded as 1 + iea n C + 0(e 2 ). For the purpose of a concise writing, we introduce 
some intermediate-step notations: 



n-1 



crt, + a. 



71- 1 



The product of the matrix elements in Eq. (2.14a) can be calculated as [23]: 



N /,Tr I iea^ — % I ,T. 



n 

n=l 



* 2 



* 2 



f N 

n 

vn=l 



(* Z J* 2 



* 2 



l* 2 



cxp 



ft 2 



AT 

E (4.-, 



n=l 



n-l 



(2.15) 



where the product of series (* z J* Zn _i) can be expressed as 



N 



n 11^.1111^. ii =cxp 



n=l 



• exp 



- Zn I 
AT 



b N+Pl N ° 2 N + l° 2 Nl h2 4>l+Pl o\ol/h 2 



2{a 



N+l 



+ <7 



N> 



2W + fS) 



n=l 



2(^+1 +<T„)^n- 



— - (c„+l + (7 n -l 



,gn+l-£n-l J,2 



1 4(a2 +l(T 2 (T 2_ i + ^2^)^ 



f N 

n 



2<J n (<T n+1 + CT„_l) 



/ 2q„cr„_i 

CT n + °n-l 



Pn-1 



£V+1-£V»-1 „2 



4n 2 



+<£)(<* 



-)]• (2-16) 



Here we introduced a virtual width <7n+i by hand, satisfying <tjv+i — <7n — &n — &n-i, m order to get the tidy sum in 
the exponential position. In the limit of N — > 00, <7n+i will approach cr/v = 07 and hence docs not effect the quantum 
dynamics. 

For the gravitational part, careful calculations outlined in the Appendix yield 



<*de- i£Q » § |*c»-. 



n=l 



TT <*C-l*C-*> cxp 

11 II*. Jill*. ,11 xp 



\n=l 



d 2 , 



ieo:„ • 37rG 
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x E (7 2 + f ) ( » (1 - d2+rf2 ) + d2+rf2 ) + sin (46„) 



2d 2 



n=l 



d n + d n _i d n + d n+1 • 



1+7 2 



((^ + f)(sin 2 (4M(l 



32 



-) + 



16 



dn + d n+1 d\ + d n _ 1 



+ 2iv n sin(86 n ) 



2dl 



■(1 



32 



dn + d n _ 1 d\ + d 2 l+1 



(2.17) 



where ' 



N 

n 



2didt 



and 



*Cnl*C„-! 



n=l 



l*fJIII* 



exp 



Cn-l 



A r 



4- ft 2 d 2 d 2 7i 2 -I- h 2 d 2 rl 2 



2(<& +1 +<&) 2(rff+^) 



' AT 

n 



■ exp 



AT 

[■EC 

71=1 



2(d 2 n+1 +d 2 n )v, 



V n — V n — l 



(d n 



+ 1 



2(d 2 l+1 +d 2 l )(d 2 l +d 2 n _ 1 ) 



h ib n 



Md 2 l+1 d 2 n d 2 n _ x + djdl^bj^^ + 2dj{d n+1 + d n ^) d ^-/^ b i 

m+i+dDK+di^) 



(2.18) 
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Now we take the limit N — > oo and substitute Jg dr for X)^Li e to S et tne functional integral formalism of the 
amplitude: 



where 



i4([* / ][* j ]) = ei(l*/l 2 -l*'l a +IOI 2 -Kir) y p a y[p0(r)][D^(T)][Dw(T)][D6(r)]( 



i(sr"+sr T ) 



JO 



dr i- 



2^2 



V4<T 2 / dr\ Ah 2 H h 2 



2<r 2 



(2.19) 



(2.20) 



d f f_\ .d_ f d 2 b 2 
Ad 2 ) +l dr\ A 



bv 



7 2 (, 2 + d -) ((sin 2 (26) (1 - 1) + I) + - sin (46) (l - ±)) 



1 + 7 2 



(„ 2 + y) ((sin 2 (46)(l - + J) + 2^ sin (86) (l - ^))]). (2.21) 



Here the "dots" over <j> and v stand for the time derivative with respect to the fictitious time t. The functional 
measures are defined on continuous paths by taking the limit of N — > oo: 



AT 



/ 2?a := li 

J N- 

/^, 1MT ,,.,„(n^:)/n 



\ f a 2o \ f a ~io 



2nh 
dv n db n 



2tt 



(2.22a) 
(2.22b) 
(2.22c) 



Ignoring the total derivatives with respect to r in Eqs. (2.20) and (2.21), we can read out the total effective Hamiltonian 
constraint as: 



^ cS = 2^ ~ inGl 



(, 2 + y) (sin 2 (26) (1 - ) + |) + - sin (46) (1 - ±) 



+ 



3^G(l + 7 2 ) 



(, 2 + f ) ( sin 2 (46) (1 - f 2 ) + A) + 2iv sin (86) (l - f ) 



16, 



. (2.23) 



Note that ^- and -tj are the square of fluctuations of and sin (26) respectively. They can be seen as quantum 

corrections to the leading term: v 2 sin 2 (26) + iv sin (46) of the Euclidean part 0e as well as the leading term: 
v 2 sin 2 (46) + 2iv sin (86) of the Lorentz part. 

A careful observation reveals that the real and imaginary parts of the leading terms can be synthesized into a Moyal 
^-product [22], i.e., 

v 2 sin 2 (26) + iv sin (46) = ve 1 9 9 ") ( sin (2b)v sin (26)) =: t; * ( sin (2b)v sin (26)) , (2.24a) 
v 2 sin 2 (46) + 2wsin(86) = v * ( sin (46)u sin (46)) . (2.24b) 

Therefore the effective Hamiltonian constraint takes the form: 



+ 



3ttG(1+7 2 ) 



v * ( sin (2b)v sin (26)) (l - ^ 
v * (sin (46)u sin (46)^ (l - + 



sin 2 (26)rf 2 , 4 , 2v 2 

16, sin 2 (46W 2 16 , 8v 2 / 

* 1 — — - — (1 I H h 4 

2 ^ d 2 J + d 2 + 



. (2.25) 



The Moyal ^-product emerges in the gravitational part of the Hamiltonian, since both 6e oc t)(sin (26)z)sin (26)) 
and Ol oc v(sin (46)£>sin (46)) are non-symmetric operators which can be regarded as a product of two self-adjoint 
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operators. Thus, in this model the coherent state functional integral also suggests the Moyal ^-product to express 
the effective Hamiltonian for the quantum system with a non- symmetric Hamiltonian operator. However, since the 
Moyal ^-product originates from the non-symmetry of the operator, one may doubt why we did not use a symmetric 
operator from the very beginning. To understand the motivation of the non-symmetric operator 9, we recall that 
the initial Hamiltonian constraint operator in LQC is actually self- adjoint in the kinematical Hilbert space [8]. To 
resolve the constraint equation and find physical states, one feasible method is to rebuild the constraint equation as a 
Klein-Gordon like equation and treat the scalar (f> as an internal time. As a result, the constrained quantum system 
was recast into a relativistic particle whose dynamics is govern by a Klein- Gordon like equation with an emergent time 
variable [9]. The price to get this Klein-Gordon like equation is that the new gravitational Hamiltonian operator 
becomes a multiplication of two self-adjoint operators, and hence it is no longer symmetric. But this does not indicate 
that one could not employ 6 in the intermediate step to find physical states. On the other hand, because the Moyal 
^-product comes from the expectation value of the multiplication of two self-adjoint operators on coherent state, this 
non-symmetric just provides a suitable arena to examine the Moyal ^-product from the path integral perspective. 

We can also take another practical way to symmetrize at the beginning. For example, one can define a symmetric 
version of Oe by 

@e : = ^(@e + ©1) « [v(s\M2b)vsiM2b)) + (sin(26)usin(26))t)], (2.26) 

and then carry out the same procedure of above coherent state functional integral. In the calculation of matrix element 
(^cJ©eI*C™-i)' we cou ld think that the operators v and sin (26)vsin (26) in 0e act on bra (\I>£ n | and ket \^c n _ 1 ) 
respectively, while sin (26)vsin (2b) and v in 0^ act on bra (\I>£ n | and ket \ i f?£ n _ 1 ) respectively. Then it is not difficult 
to see that the imaginary parts generated by ©e and 0g cancel each other. Hence for the symmetric Hamiltonian 
operator corresponding to 0, we can finally get the following effective Hamiltonian constraint: 

U := - ^ + 37rG (f + ^2) ( sin2 ( 25 )(! " ( 16 + 12 7V - (1 + 7 2 )(1 - 16e 2 ) sin 2 (26)) + 2e 2 ) , (2.27) 
which takes the same form as (2.25) but without ^-product while e = 1/d denotes the quantum fluctuation of sin 6. 



III. EFFECTIVE DYNAMICS 



Using the effective Hamiltonian constraint Jif e s which contains Moyal ^-product, one may investigate the corre- 
sponding dynamics by defining the evolution equation as: 



f(v,b) :=-(/*J^eff-^eff*/), 

for any dynamical quantity f(v, b). Especially, the evolution of basic variables can be obtained as: 



v = — - 



127rG 7 2 



h 



v * (v sin (26) cos(26)(l-4e 2 )) + 



2 , sin (26) cos (26) (1 - 4e 2 ) 



2e 2 



v 2 ( , 1 \ . 2 sin 2 (26)(l-4e 2 ) 
~2~V + 2?' Sm {2b) ~ 



8e 4 



d h e 2 



+ 



127rG(l+7 2 ) 
Ah 



2v * (v sin (46) cos (46) (1 - 16e 2 )) + 



2 , sin (46) cos (46) (1 - 16e 2 ) 



sin 2 (46)(1 - 16e 2 



6 = 



+ (V-(, 2 + ^)4sin 2 (46)- g£4 
' T( ; [ (2v(l - 4e 2 ) sin (26)) * sin (26) + Ave 2 



d b e 2 



_ ( ,sin 2 (26) (i _ 4g2 , 



2e 4 



+ ^L)4sin 2 (26) -2v 2 y v e 2 



3^G(l + 7 2 ) 



Ah 



(2v(l - 16e 2 ) sin (46)) * sin (46) + 16ve 
sin 2 (46) 



2e 4 



1 - 16e 2 ) + (v 2 + 16 sin 2 (26) - 8v 2 jd v e 



(3.1) 



(3.2a) 



(3.2b) 
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where dbS 2 = d(e 2 )/db and d v = d/dv. However, there seems no way to understand Eqs. (3.2a) and (3.2b) directly 
as effective classical equations because of the ^-product therein. Instead, we could use the effective Hamiltonian 
constraint (2.27) without ^-product to explore the effective dynamics. Using the conventional Poisson bracket, we can 
derive a modified Friedmann equation from the effective Hamiltonian (2.27) as: 



H 2 = 



a\ 2 8irGp c 



2v 2 e 



(sin (26) cos (26) (1 - (16 + 12 7 2 )e 2 - 2(1 + 7 2 )(1 - 16e 2 ) sin 2 (26)) 
- ( sin 2 (26) (4(1 + 7 2 ) cos 2 (26) - 7 2 ) - ±) d b e 2 ) 
- ( sin 2 26(1 - (16 + 12 7 2 )e 2 - (1 + 7 2 )(1 - 16e 2 ) sin 2 (26)) + 2e 2 J (3.3) 

where p c = s^^Wy 5 * s a cons tant. To annihilate sin (26) and cos (26) in Eq. (3.3), we use the constraint equation 
(2.27) to get 



and 



. 2 1 - (16 + 12 7 2 ) £ 2 - ^(1 - (16 + U^e 2 ) 2 - 4(1 + 7 2)(1 - 16 £ ^ 
Sm {2b) = 2(1+ 7 2 )(1-16 £ ) 2 ' (3 ' 4) 

x ^^-2 £ 2 , (3.5) 

J Pc 

where p — is the density of matter, K = 1 + 2 ^ 2 t , and J = 1 + 2 J> ei . However, Eq. (3.3) looks problematic since 

it depends on the volume v of the chosen fiducial cell. This originates from the fact that we have to use the coherent 
states peaked on the phase points (v, 6) in the path integral. In the final picture we have to remove the infrared 
regulator by letting the cell occupy full spatial manifold. In this limit, the irrelevant correction terms proportional to 
l/(ve) 2 could be neglected, while the relevant terms proportional to e 2 would be kept, since e was understood as the 
fluctuation of sin 6 which does not depend on the fiducial cell. We finally get 



h2 = 8nGp c 



( A - V / A 2 -4/3 X ) (2/3-A +v /A2-4/3 x ) (a 2 - 4/3 X ) 

( 1 A - y/A 2 - 4/3 X , 4(1 + 7 2 )(2/3 - A + y^A 2 - 4/3 X ) 2 ,x 2 ^ 
+ V2 20 ( 20 7 ) r b£ 



(3.6) 



where the positive and negative signs correspond to the expanding and contracting universe respectively. Here we 
use notations: A = 1 — (16 + 12 7 2 )e 2 and = (1 + 7 2 )(1 — 16e 2 ) for a concise writing. Note that Eq. (3.6) implies 
significant departure from classical GR. For simplicity, we first consider the case dbS 2 = and see whether the bounce 
or re-collapse determined by H = could occur. Then it is obvious that, for a contracting universe, the so-called 
quantum bounce of LQC will occur when 

X=%, (3.7) 

which means that p increases to pboun ~ 4(^^,2) if £ 2 is neglected. On the other hand, for an expanding universe, 
a re-collapse would occer when \ — 0, or equivalently p decreases to p co \\ « 2e 2 p c , which coincide with the result 
in canonical theory [18], where d^e 2 was assumed as higher order term and hence neglected. As pointed out in 
Refs.[17, 18], the inferred re-collapse is almost in all probability as viewed from the parameter space characterizing 
the quantum fluctuation e. Intuitively, as the universe expands unboundedly, the matter density would become 
so tiny that its effect could be comparable to that of quantum fluctuations of the space-time geometry. Then the 
Hamiltonian constraint may force the universe to contract back. It should be noted that the effective equation and 
hence the inferred effect of re-collapse are only valid with the coherent states of Gaussian type. Whether there is a 
similar result for other semiclassical states is still an interesting open issue. For example, one may consider the affinc 
coherent states in the affine quantum gravity approach developed by Klauder [26, 27]. 

In the case when d^e 2 could not be neglected, the bounce would also be approached for a contracting universe. 
Because d^e 2 would not be bigger than the order of 0(e 2 ), x could be infinitely close to the result in Eq. (3.7) and 
lead to H = 0. However, for an expanding universe, both of the two terms in the bracket of the right hand side of 
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Eq. (3.6) are non-negative in large scale. As a result, the Hubble parameter would always keep non-zero unless dbS 2 
approaches asymptotically. Therefore, the inferred re-collapse might occur only if dt,s 2 approaches asymptotically. 

If we neglect all the higher-order quantum corrections: -^2, e 2 and dbS 2 , Eq. (3.5) would be simplified to \ = y, 
and hence a first-order modified Friedmann equation could be obtained from Eq. (3.6) as: 



H 2 = 



8irGp 



_ 7 2 + 4(1 + 7 >/pc l 2 Pc 

1 +7 2 2(l+ 7 2 ) 2 p 



1 - 



4(1+7> S 



1-4/1- 



4(l + 7 2 )p 



(3.8) 



Note that this first-order modified Friedmann equation is different from that in Refs.[8, 17]. But it coincides with 
the modified Friedmann equation in Ref.[18]. Hence Eq. (3.8) still contains the particular information of alternative 
dynamics. It is easy to see that if the matter density increase to p = ^ppy, Hubble parameter would be zero and 
the bounce could occur for a contracting universe. On the other hand, in the classical regime of large scale, we have 
X <C 1 for p <C p c and hence Eq. (3.8) reduces to the standard classical Friedmann equation: H 2 = 8nGp/3. 



IV. SUMMARY 



Since there are quantization ambiguities in constructing the Hamiltonian constraint operator in LQC, it is crucial 
to check whether the key features of LQC, such as the quantum bounce and effective scenario, are robust against 
the ambiguities. Moreover, since LQC serves as a simple arena to test ideas and constructions induced in the full 
LQC, it is important to implement those treatments from the full theory to LQC as more as possible. Unlike the 
usual treatment in spatially flat and homogeneous models, the Lorentz term has to be quantized in a form quite 
different from the Euclidean one in full LQC. For the above purpose, this kind of quantization procedure which kept 
the distinction of the Lorentz and Euclidean terms was proposed as alternative dynamics for LQC [18]. It was shown 
in the resulted canonical effective theory that the classical big bang is again replaced by a quantum bounce and it is 
possible for the expanding universe to re-collapse due to the quantum gravity effect by certain assumption. Hence it 
is desirable to study such kind of predictions from different perspective. Meanwhile, it is also desirable to study the 
Moyal ^-product by coherent state functional integral approach within LQC models. 

To carry out the above ideas, the present paper is devoted to study the coherent state functional integral in spa- 
tially flat isotropic FRW model coupled with a massless scalar field in the alternative dynamics framework of LQC. 
The main results can be summarized as follows. By the well-established canonical theory, the coherent state func- 
tional integral for LQC with alternative dynamics has been formulated by group averaging. For the non-symmetric 
gravitational Hamiltonian constraint operator, the Moyal ^-product emerges naturally in the resulted effective Hamil- 
tonian with higher-order quantum corrections. For the corresponding symmetrized Hamiltonian operator, the effective 
Hamiltonian and modified Friedmann equation are also derived from the coherent state functional integral approach. 
It turns out that the quantum bounce resolution of big bang singularity can also be obtained by the path integral 
representation. On the other hand, if higher order corrections are included, there is a possibility for the re-collapse 
of an expanding universe due to the quantum gravity effect, which coincides with the result obtained in the canon- 
ical formalism. Moreover, the first-order modified Friedmann equation still contains the particular information of 
alternative dynamics and hence admits the possible phenomcnological distinction between the different proposals of 
quantum dynamics. The alternative modified Friedmann equation (3.6) or (3.8) sets up a new arena for studying 
phenomenological issues of LQC. 



Acknowledgement 

We thank Abhay Ashtekar and Peng Xu for helpful suggestion and discussion. This work is supported by NSFC 
(No. 10975017) and the Fundamental Research Funds for the Central Universities. 



9 



Appendix A: Calculation of the functional integral 



We give some details on the calculation of the Lorentz part of the matrix element of exponentiated gravitational 
Hamiltonian operator: ( l J / cJ e_i£aeL I^C^-i) ■ The order of 0(e) of this matrix element is 

(* c „ |-*eae L |*c»-i) = ^ ^^V 7 ^ S [ k( * k + ( fc )*C.-i (fc + 8) - 2fc 2 *^ (fc)*C„-i (fc) 



3ttG(1 + 7 



+fc(fe-4)*J B (fe)* CB _ 1 (fc-8) 

| (j (-^n.n-l — -^n.n-l + ^n,n-l) ■ 

Now we need to deal with the three terms L^ n _ 1: LP n n -i,L~ n _ l separately. First, we get 

(fc-„„) 2 (fc+s-^„-i) 2 



8{v„ - v n -i) 



2di 2d 



: g — »6n(fe— Wn)+»i>n-l(fe+8— Vn-l) 



k 

— exp 



32 j- (v n - Vn-l) 2 r 

+ i&K - —, , n ; , 9 — - + «o„(w„ - V n -l 



4 + 4-i € + " 2(^ + 4^) 

J> 2 + 4fc) exp (-(fc - tJ+) 2 /4 - i(&n - b n -i)(k - TJ+)) 



(Al) 



(A2) 



where v 



. To do the summation in the above equation, we first have to rewrite k 2 + Ak as a 



function of fc — v,. 



k 2 +4k = (k-v+y + (2v++4)(k-v+)+v+ +4^ 



(k-v+) 2 + 2(v n - 



64 - 2dj_ 1 
4 + 4-1 



(k - v+) + v n 2 - v n 



2 _ 12#-4(£_i , 32d 2 (d 2 -rf 2 _!) 



+ 



4 + 4-i (4 + 4-i) 



2 \2 ' 



where u„ = d " 1 • Do the summation term by term, we get 



_ 3Z _ I 



8< 4 + \ 

d 2 +4-i 2 +^.»-iJ' 



(A3) 



where P^n-i denotes a polynomial of u„ — u n _i, b n — 6 n _i and d„ — rf„_i without the zeroth order term. Here we 

have expanded the factor exp m Eq. (A2) as 1 — 4 ^ , 2 "~J 2 "~ 1 ' > H • Except for the leading term 1, all 

the other terms can be conflated with P^n-i- Under the continuous limit N — > 00, this P^ n _ 1 does not contribute 
to the effective action of gravity. It is easy to calculate lP n n _ 1 and L~ n _ 1 as follows: 



L° n , n -i = 2<*c»i*C„-x> ((«™) 2 + f +P°, n -i 



Kn-i = (*C„l*C.- 1 )e" 3ST ^ I7 e- i85 « (v n f +v n 



Ml 



dl 



d 2 n + 4-1 2 



(A4) 
(A5) 



Taking the expansion e 
get the combination 



( <+k-i) = 1 



d 2 + 8 d -2 hO(y) and neglecting the higher order terms than (3?), we can 



r+ -7"° _i_ r- 

-^n.n-l n n,n-l ' ^n,ra-l 



= -4<* Cb |* Cb _ 1 ) [((« n ) 2 + f ) (sin 2 (4h„)(l 



-) + 



16 



4 + 4V 4 + 4-i 
+ 2,sin(86,0,„(l-^^)^^) + P„T-i 
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and hence the matrix element (\&f n |e * eaeL |^'^ n _ 1 ) is 

. 37rG(l+ 7 2 )/ /fc ,2 , iftw, 2 



(*<„l*C„-i>«»P 



(((5o)! + f)(si », (4tn)(1 __J^ ) + _^_ ) 

32 2cZ 2 
+2. sin (86„K (1 - -^—^—) + ^ 



(A6) 



Similarly, we can get the Euclidean part (^cje " a6>E I^Cn-i) as: 



(^Je-^l^n-x) = (^J^Cn-J exp [^„3^G 7 2 (((tJ„) 2 + ^)(sin 2 (26„)(l - rf2 + " rf2 ) + rf2 + 4 rf2 ) 



4 + 4V 4 + 4-1' 

+ i ™<«»)"-( 1 - dfnrr'rtT + P " B "-')1- <A7) 



Combining Eqs. (A6) and (A7), we can get Eq. (2.17). 
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